Collective charge excitations along cell membranes 
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A significant part of the thin layers of counter-ions adjacent to the exterior and interior surfaces of 
a cell membrane form quasi-two-dimensional (2D) layers of mobile charge. Collective charge density 
oscillations, known as plasmon modes, in these 2D charged systems of counter-ions are predicted in 
the present paper. This is based on a calculation of the self-consistent response of this system to 
a fast electric field fluctuation. The possibility that the membrane channels might be using these 
excitations to carry out fast communication is suggested and experiments are proposed to reveal 
the existence of such excitations. 

PACS numbers: 87.16.Dg,87.14.Cc,82.39.Wj,82.45.Mp 



Collective charge excitations in condensed matter, the 
so-called plasmon modes, play a fundamental role in un- 
derstanding the properties of a metal 1 . These are exci- 
tations of the electron gas as a whole and are respon- 
sible for the electromagnetic response of the metal for 
frequencies above and below the classical plasma oscil- 
lation frequency, and for our understanding of other ef- 
fects such as screening. The plasmon is the quantum 
of such collective charge oscillations which may exist in 
any charged system and in fact they were first discov- 
ered in a classical plasma 2 *. In quasi-two-dimensional sys- 
tems such as surfaces of metals^^^, or interfaces of 
metals with insulators, plasmons have been investigated 
both theoretically^iSiLS and experimentally^ for several 
decades. 

Plasmon modes have been found^ in a system of elec- 
trons on the surface of liquid helium. The surface of 
a lipid bilayer with the surrounding counter-ion system, 
which we consider in this paper, is analogous to that of 
electrons on helium surface as will be made clear be- 
low. The surface of a cell membrane in aqueous en- 
vironment becomes negatively charged^. It is well- 
knowni2iiiii2ii&±i that a diffuse double layer of counter- 
ions from the solution screens the negative charge, thus, a 
thin two-dimensional (2D) layer of mobile cations is accu- 
mulated adjacent to the extra-cellular and intra-cellular 
membrane surfaces. In addition, because of different 
permeabilities of the cell membrane channels to various 
cations, the inside of the cell is kept at a negative poten- 
tial relative to the outside. In this paper it is pointed out 
that the thin layers of counter-ions adjacent to the lipid 
bilayer are characterized by 2D longitudinal plasmon ex- 
citations. These modes are excited in response to a fast 
perturbation such as the opening of an ion channel. The 
gating mechanism of ion channels along the cell mem- 
brane, that are believed to open using voltage sensing 
gates, is an open problem in cell biologyi£*i&. The collec- 
tive charge density wave which is theoretically found in 
this paper can play the role of a very fast signal carrier 
for such ion channel communication. Therefore the col- 
lective modes investigated here can play a fundamental 
role in cell biology. 



The negatively charged exterior and interior surfaces 
of biological membranes of cellaii as shown in Fig. 
where the density of charged phospholipid heads is higher 
on the exterior membrane surface. Diffuse double lay- 
ers of counter-ions, such as H + , Na + , K + , Ca ++ , 
form on each membrane side. Due to the different mem- 
brane permeability for these cations, there is a voltage 
— V (of the order of — lOOmV) between the cell inte- 
rior relative to the exterior. The solid line in Fig. ^ is 
the positive free counter-ion distribution near the mem- 
brane surfaces^. In thermal equilibrium the counter- 
ion charge density varies as we move away from the 
membrane surface as shown, within the so-called Gouy- 
Chapman(GC) length X GC = kBTei/2iraiqi for each 
membrane side. Here q and qi are the dielectric con- 
stants and the counter-ion charges of the outside and 
inside of the cell. The surface densities ai — qin® of 
both membrane sides are experimentally accessible and 
they can be as high as n° ~ 1Q~ 2 A~ 2 or even higher. 
This counter-ion surface charge distribution is controlled 
by the bulk cation concentrations of the intra-cellular 
and extra-cellular fluids^. In addition, depending on the 
temperature, cation concentration, and cation charge, a 
fraction of these ions can be adsorbed on the membrane 
surface. The adsorbed ion density versus bulk aqueous 
ion concentration is determined from the Stern adsorp- 
tion isotherm theory 14 and its extensions 1819 . The ap- 
proximate value of Xgc is about A A, if we use T = 300K, 
e = 80 and n° — 10 A . For charge fluctuations prop- 
agating with wave vectors k parallel to the surfaces and 
magnitude k << 1/Xgc (which is the case of our inter- 
est), these layers can be regarded as two-dimensional. 

In this paper, the collective charge response of the mo- 
bile counter-ion charge of the interior and exterior surface 
layers to an external electric held perturbation is studied. 
It is shown that there are longitudinal charge collective 
modes with wave vectors parallel to the surface and the 
two layers are coupled and fluctuate together as a whole 
in the limit where the wavelength is much longer than 
the membrane thickness d. 

The counter-ions of the diffuse double layer adjacent 
to the extra-cellular membrane surface are character- 
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ized by mass mi charge q\ = Z\e [Z\ the ionic va- 
lence and e the electron charge) and instantaneous sur- 
face number density ni(r,i) (r = (x,y)) and those ad- 
jacent to the intra-cellular membrane surface by mass 
7712, charge q 2 = Z 2 e and instantaneous surface number 
density 772 (r,i). We are interested in studying devia- 
tions 5ni(r,t) from the equilibrium surface density rq, 
where 5rii(r, t) — 77.; (r, t) — nf. Following the well-known 
approach 1 to this class of problems where there is the 
possibility of collective response of the charge, a pertur- 
bation </> ext (r, t), due to an external charge density fluc- 
tuation p ext (r, z,t) and current fluctuations 3 ex t(r, z,t), 
is introduced. These physical quantities are related as 
follows 



dtpext(r,z,t) +V 



= -4irp ext (r, z.t) (1) 
3 ex t(r,z,t) = 0, (2) 



where the z direction is taken perpendicular to the mem- 
brane surface. This perturbation may be caused by the 
opening of a sodium or potassium membrane channel or 
by a real external probe. While these channels are part 
of the cell membrane, they are not part of the degrees 
of freedom considered here, namely we only consider the 
system of the accumulated charge on both sides of the 
cell membrane. 
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FIG. 1: The lipid bilayer cell membrane of width d and its 
various features are explained in the text. 



We would like to investigate the response of the sys- 
tem to such an external perturbation namely to solve the 



Poisson equation: 

V 2 ^(r,z,i) 



4?r 



p(r,z,t), 



(3) 



where i = 1,2 and denotes the cell exterior and inte- 
rior respectively. Here 0i(r, z,t) is the total field in the 
regions i = 1, 2 and the total charge density is given by 

p(r, z, t) = po(z) + p ext {v, z, t) + p md (r, z, t) (4) 

where po(z) is the equilibrium charge density and 
Pmd(r, z,t) is the induced charge density. Within our 
2D approximation, the induced density is written as 

Pind(r,t) = qi5n 1 (r,t)5{z) q 2 5n 2 {Y,t)5{z + d). {h) 

For both layers we require the validity of the linearized 
continuity equation 



dtSm{r,t) +77°V • Vi (r,t) = 0, 



7 = 1,2 (6) 



where Vi(r, t) is the velocity of the counter- ions along the 
i th layer. The equation of motion for a single counter-ion 
maybe written as follows: 



,; (d t Vi + — 



- qiV<f>i(T,z = Zi,t) 

> 

-VSrii(r,t) 



rriisf. 



(7) 



where Zi — CP,<i + for the exterior and interior layer re- 
spectively and r is the time between collisions, while 
Si = \J ^UbT ' jrni is the velocity of density fluctuations 
in the charged subsystems and the value of 7 = 2 will be 
used for a 2D ideal gas. The first term is the force on 
the counter-ion due to the local electric field while the 
second term is the force due to the local compression of 
the counter-ion 2D gas. In order to simultaneously solve 
Eqs. (|3lt)l7fl . we consider the Fourier transform (denoted 
by a tilde) of all the functions of the 2D position vector 
r and time t. Eq.© takes the form: 



An 



(d 2 z -k 2 )4> t (k, Zl u;) = -— p(k,z,Lj) (8) 



where k is the wave vector parallel to the membrane sur- 
face and ui is the frequency. Eq. (JHJ has the following 
solution: 



2tt 



[qxSnx (k, u)e~ k ^ + q 2 6n 2 (^ w)e~ fc|2 ~ d| 
' dz'p'(z')e- k \ z - z '\]. (9) 



where p'(z) = Pq(z) 



Pext 



(k, z, u>). Substituting J5J in 



the Fourier transform of Eqs. H6I7|) an d by eliminating 
the velocity variables we find 



oj 2 {k) (?77i(k, uj) — Pi(k)Snj(k : uj) 
7i(k,w), (10) 



where i,j — 1,2 or 2, 1. We have used the definitions 



A(k) 



^ I dz'p'(z>)e-V 



Z 1 Z 2 Q 2 i e- kd , 
Z 2 uj 2 + s 2 k 2 , 
2ire 2 n% 



(11) 

(12) 
(13) 

(14) 



The system of Eqs. (|1(J[) has a solution for the responses 
<5n;(k,a->) which attain their maximum values for the fol- 
lowing frequencies: 



1 N - 

27/ 



n±(fe) = \/<4(fc) + (£) +^ 



(15) 



c4(fc) 



w 2 (fc) 



± 



u> 2 {k) - u 2 {k) 



f3 1 {k)(3 2 (k). (16) 



The frequencies uji(k) are the frequencies of plasmon os- 
cillations of each layer when their mutual coupling (terms 
(3i(k)) is small, i.e., when kd >> 1. 

The frequencies are shown in Fig. [5] for the case where 
all the parameters of both layers are taken to be the same. 
The frequencies are given in units of 



w 



27re 2 n° 
mde 



(17) 



and the wave- number in units of 1/d. In terms of these 
variables the frequency of each independent layer can be 
expressed as 



! (fc) 



= kd(l + kd ) 



d Q = k B Te/{2Tr 2 e 2 n°) ~ 5A 



(18) 
(19) 



and the second term has a small contribution for the 
range of k in the graph because d ~ 30A The frequencies 
of each independent layer are shown by the middle curve. 
The lower frequency corresponds to in-phase charge os- 
cillations of the two layers and the higher frequency to 
out-of-phase charge oscillations. It is straightforward to 
show that in the long-wavelength limit (kd << 1) 



LU-(k — > 0) = C p k 



u + (k ->0)= ^2s 2 k/d+{s 2 - s 2 )k 2 



(20) 
(21) 



and in the particular case where the parameters char- 
acterizing the layers are equal, we have c p — \A s o + s i 
with so = wod. Taking realistic values of the parameters 
we find that for H + counter-ions the frequency-scale is 
wo/27r ~ 10 u Hz and the speed is c p ~ 3000to/s. For 
heavier counter-ions these scales will be smaller as both 
scale as 1/y/m. 

The group and phase velocity of the higher frequency 
mode and of the independent layer oscillations in the 




FIG. 2: Eigen- frequencies of the membrane collective charge 
oscillations when all characteristic parameters of the carriers 
of both layers are the same. The curve in the middle denoted 
by lji is the plasmon frequency of each layer in the limit where 
the coupling between the inner and outer layers is neglected. 
The frequencies denoted by u- (w+) correspond to in-phase 
(out-of-phase) oscillations of both layers together. When the 
parameters of the two layers are different, the corresponding 
plot looks similar, where instead of having one Wi curve in 
the middle we have two, one for Wi and a second for ob- 



long-wavelength limit diverge as u(k — > 0) = s^/ykd. 
Therefore in the long-wavelength limit these velocities 
are higher than the counter-ion thermal velocities and 
therefore the collision frequency is lower than the collec- 
tive mode frequency. Therefore the Landau damping rate 
is small in the long- wavelength limit. In order to obtain 
a more accurate expression of the plasmon frequencies 
instead of the bare mass and surface density, one should 
use the effective mass of the hydrated ion and the effec- 
tive nearly-free counter-ion charge density. If there are 
binding sites on the membrane surface this will reduce 
the effective free surface charge density. Therefore more 
realistically we should expect frequencies below the above 
mentioned values. In addition, the surface corrugations 
of the electrostatic interaction of the counter-ions with 
the phospholipid surface introduce a non-zero frequency 
gap Qq at zero plasmon momentum. This fio is the fre- 
quency of local charge oscillation near a local minimum 
of the interaction potential between the counter-ion and 
the membrane surface. When the surface charge density 
is low, this interaction may lead to the formation of a 
Wigner crystal where the counter-ions will order in a tri- 
angular or other lattice structures^. These structures of 
ordered charged layers have been observed in electronic 
bilayer systems, consisting of two quasi-two-dimensional 
layers of electron or hole liquids^. Such a Wigner crystal 
will be characterized by plasmon modes with frequency 
given by the same form as Eqs. I|15lltj|) with the masses 
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Lipid bilayerfilm'XS;Ssi 

FIG. 3: An idealized experimental set up to control the cation 
concentration on both membrane sides. 



replaced by the ion effective mass within the Wigner crys- 
tal. 

These collective modes can be investigated experimen- 
tally by studying the electromagnetic absorption spec- 
trum in the frequency region of 10 10 — lO 11 ^^ (de- 
pending on the wavelength of the excited plasmon and 
the counter-ion parameters and distribution) of artificial 
phospholipid bilayer membranes which separate aqueous 
solutions of different ionic concentrations (see Fig. 
with high counter-ion concentrations^ <Ji near the ex- 



terior and interior surfaces. Up to 10mm 2 size mem- 
brane can be formed in the hole of the wall separating 
the two aqueous compartments shown in Fig. by ap- 
propriately mixing the phospholipids as was first discov- 
ered by Mueller et alw An appropriate voltage can 
be applied between the two aqueous compartments of 
different cation solutions separated by the membrane. 
Plasmon modes of the two thin counter-ion layers ad- 
jacent to the membrane can be excited and probed by 
a source/spectrometer system. With the experimental 
techniques and source/spectrometer systems available to- 
day, it is possible to study the AC electromagnetic re- 
sponse of very small samples over a range of frequen- 
cies and to cover the above frequency domain with one 
or a combination of source/spectrometer systems. The 
plasmon frequencies should scale as yfa with respect to 
surface concentrations (which can be estimated from the 
bulk concentrations) and as 1/ y/m* with respect to var- 
ious counter-ion effective masses. 

These 2D plasmon oscillations can be used by the var- 
ious parts of the cell membrane to carry out commu- 
nication. For example, these modes are excited in re- 
sponse to a cell membrane pore opening which creates 
a local charge distribution fluctuation. The phenomeno- 
logical studies of the behavior of the gates of the voltage- 
sensitive ion channels in cell membranes indicate a high 
sensitivity^ of the voltage sensor to small changes of the 
membrane electric field. This implies that the 2D plasma 
oscillations could play a role in triggering the sensor to 
open the gate. Before, however, one proceeds further to 
study such possibilities it is important to carry out an 
experimental study to verify the existence and the prop- 
erties of these plasmon modes. 
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